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Abstract. Residual interactions in many particle systems lead to strong cor- 
relations. A multitude of spectacular phenomenae in many particle systems are 
connected to correlation effects in such systems, e.g. pairing, superconductiv- 
ity, superfluidity, Bose-Einstein condensation etc. Here we focus on few-body 
bound states in a many-body surrounding. 

1 Introduction 

Correlated many particle systems such as, e.g. nuclear matter or strongly cou- 
pled plasmas, have a complicated dynamical behavior. Only few areas of the 
density-temperature phase diagram used to characterize the state of the sys- 
tem in thermal equilibrium may be described in the approximation of nonin- 
teracting quasiparticlcs. The dynamics of the quasiparticle is determined by 
the mean field of the other particles and therefore some prominent features like 
self-energy corrections and Pauli blocking (Bose enhancement) are sufficiently 
taken into account |^ . However due to sizable residual interactions many inter- 
esting and exciting phenomena, such as clustering, formation of condensates, 
and phase transitions occur. The description of the dynamics of such systems 
is particularly difficult since the quasiparticle approach reaches its limits. 

A large number of these phenomena such as, e.g. two-particle bound state 
formation, can already be accounted for by explicitly introducing two-body 
correlations into the formalism. This may be achieved in the frame work of 
the Green function method 0| and leads to effective two-body equations that 
include medium effects in a consistent way. 
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Figure 1. Density temperature pliase 
diagram of nuclear matter, no denotes 
nuclear matter density no = 0.17 fm^''. 



Figure 2. Nuclear matter phase dia- 
gram. The the amount of correlated den- 
sity is given in per cent of the total den- 
sity Tc denotes the critical tempera- 
ture. 



The density temperature planes of nuclear matter are shown in Figs. |l| and 
^. The phase diagram of nuclear matter turns out to be very rich. In particular 
the superfluid phase reflecting strong pairing is relevant for the structure of 
neutron stars Q. At lower densities bound states occur. This part of the phase 
diagram may be accessed in the laboratory through heavy ion collisions at 
intermediate energies. The description of these phenomena in the framework 
spelled out here constitues a wide and exciting area for an application of few- 
body methods. 

Besides the pairing of nucleons, correlations with more particles are also 
important. Based on the suitably modified AGS (Alt-Grassberger-Sandhas j^) 
formalism the reaction cross section for break-up {Nd — > NNN), relevant e^. 
for the formation of deuterons in a heavy-ion reaction, has been given in 
As a first application the life time of deuteron fluctuations in nuclear matter 
was calculated in |^]. In we have given an exploratory calculation of the 
four-body condensate in strongly coupled fermionic systems. 



2 Theory 

The formalism to derive few-body Green functions within the Dyson equation 
approach or cluster mean field approximation at finite temperatures and den- 
sities has been given elsewhere Here we give some of the basic results and 
focus on bound states. 

Let the Hamiltonian of the system be given by 

H{t) = J2HoiM{t)ai{t)+ I ^ 1^2(12, l'2')aIW4(t)a2Kt)ai'W (1) 

1 121'2' 

where ai{t) etc. denotes the Heisenberg operator of the particle with quan- 
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tum numbers si, fci, etc. for spin, momentum etc. In the equal time for- 
malism we may introduce creation and annihilation operators Aa(t), where 
Aa — fli, a|a2, 010203, etc. for one-particle, hole-particle, three-particle opera- 
tors. We describe the propagation of these clusters via the Green's function at 
T = 0. However, generalization to T 7^ is straight forward using the Matsub- 
ara technique [Q . The definition of the Green's function reads 

Gl/ = -t{TA^{t)Al{t')), (2) 

where averaging is due to the density operator pq. For an open system in 
thermodynamical equilibrium the extremum condition for the entropy leads to 
the following expression for the density operator of a grand canonical quantum 
ensemble, 

The temperature (1//3 = T) and the chemical potential fi are the corresponding 
Lagrange parameters. The Dyson equation is 

^dtGi;' = 5{t-t'){[A^,Al]^)+Y, fdtM'-'gl-;'. (4) 

We use Afap = {[Aa, A^^]±) in the following, where the upper (lower) sign reflect 
Fermi (Bose) type operators. The mass matrix is given by Q 

Miz; = Sjt - t)MU,o^ + (5) 

cluster — mean — field retardation 

iM'MU,o = ([K,i7](t),4(t)]) (6) 
{M'-''M)c.pArr = -i{T[Aa,H]{t),[H,Al]{t')),,,. (7) 

For Aa = oi, i.e. the one-particle case, neglecting the term TWiir (represent- 
ing all retarded irreducible diagrams) leads to the self-energy shift in mean 
field (e.g. Hartree-Fock) approximation. For a given potential this constitutes 
a selfconsistent problem to determine /i and /3, to be solved by iteration 

/^"/(^^)- e^(-.-.) + l ' ^^ = |^+E^^(12)/.. (8) 

For r — > this function approaches a step function at the respective Fermi 
momentum. 



3 Two-body correlations, bound states and Cooper pairing 

Nuclear matter cannot be considered as a system of independent quasi par- 
ticles but for a large part it is correlated and, e.g. in the superfluid phase as 
demonstrated in Fig. M. The basis to treat correlated densities is provided 



4 

EJMeV] 




10" _j 10" 

n [fm ' ] 



Figure 3. The deuteron binding energy 
as a function of the nuclear density for 
T = 10 MeV. P denotes the relative mo- 
mentum between the deuteron and the 
medium. 
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Figure 4. Deuteron life time as a func- 
tion of the deuteron momentum at finite 
temperature (T — 10 MeV) and density 
n = no/25 |o[. 



by a generalization of the Beth-Uhlenbeck approach The nuclear density 
n = n(fi, T) is given by 

= '^ftcc + ricorr, «corr = 2n2 + Su^ + . . . , (9) 

where ri2,3 denotes the two, three-particle correlations, present as 
bound/scattering states. In first iteration these correlations may be treated 
on the basis of residual interactions between the quasi particles. The exact two 
particle equations to be solved are known as Bethe-Goldstone or Feynman- 
Galitski equations depending on some details . In ladder approximation the 
two-body equation reads 

{E - £1 ~ 62) ME) = (/1/2 - /1/2) V2 ME), (10) 
where / = 1 — /. The normalization is given by 

E ^12ihf2 - hh)^12 = (11) 
12 

where the dual wave functions are defined by 'F" = (1 — fi — /2)'f'''- Intro- 
ducing the two-body t-matrix in a standard fashion both bound and scattering 
states have been solved . In first approximation V2 of Eq. (|l^) corresponds 
to the bare force. However, at higher densities the force may be substantially 
screened. 

As an example for an in-medium bound state the deuteron binding energy 
is shown in Fig. |^. The Mott density is defined through the condition Ed = 0. 
Note that the Fermi functions /i for particle 1 etc. depend on the relative 
momentum P between the deuteron and the medium. 

The deuteron in the medium also has a finite life time. To second order 
Born approximation the life time can be extracted from the square of the 
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Figure 5. Decay of deuteron in nuclear matter induced by surrounding nucleons. 



graphs of Fig. |5|. One sees that the width Jj is given by the incoherent decay 
of the deuteron into 3p-lh states. At zero temperature one evaluates the 3p-lh 
density of states 



53p-ih = TriTr2Tr3Tr4{0(i/i - eF)e{H2 ~ eF)e{H3 - eF)9{eF - H4)} 



Since the normalization of the deuteron wave function brings a factor [E — 
2ei?)~^ we have oc [E — 2£fY- At finite temperatures T ^ ef we get as in 
Fermi liquid theory for zero sound a behavior for F^. In a more elaborated 
theory one has to consider rescattering of the incoming deuteron with the 
nucleon(s) of the medium. Summing up all the respective diagrams leads to 
a generic three-body in-medium scattering matrix, that has been derived in ^, 

As a consequence the reaction rate depends strongly on the surrounding 
medium. The numerical result for the life-time of the deuteron is shown in 
Fig. ^, where we in addition compare the use of the isolated to the in-medium 
break-up cross section at T = 10 MeV and a nuclear density of n = 0.007 
fm~^ ~ no/25 [Q, |lO). The medium dependent elementary cross section leads 
to shorter fluctuation times. As expected the difference is vanishing for larger 
momenta. In turn this effect leads to faster chemical relaxation times that is 
e.g. relevant for the notion of freeze out in heavy ion collisions. 

The appearance of a sharp pole at i? = 2/i in the respective Green function 
related to Eq. (|lO|) corresponds to the Thouless criterion for the onset of super- 
fluidity for np Cooper pairs. Depending on the channels considered, it allows us 
to determine the critical temperatures F^ or F^ for the respective isospin sin- 
glet and triplet channels, as a function of the chemical potential. This is shown 
in Fig. H and later in Fig. || for Ff. Using the canonical relation Akp = 1.76Tc 
one realizes that the np gap has a maximum value of ~ 10 MeV. This is more 
than 3 times larger than the nn gap! Such high values may have an influence 
in expanding nuclear matter as produced in central heavy ion collisions. 

An interesting extension of Eq. (|l^) is to write this equatioir down for the 
case where we have a condensate of pairs, i.e. a superfluid state. For pairs at 
rest one has 
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Figure 6. The deuteron (Cooper pair) wave function. Large component: s-wave; 
small component: d-wave. The Paris force was used in the gap equation and the 
single particle energies are from Brueckner Hartree-Fock. (From 



where Vpk is the two-body interaction in momentum space and the occupation 
numbers Up (rip — 1 ~ Up) are given by 



n„ = — 
^ 2 



1 - 



Ep- H 



(1 - 2/p 



(15) 



with Ap = Vpk^k the usual gap function. The chemical potential /i and 
the density are related by N/V — 2J2p''^p- We see that for a superfluid state 
Eq. (|l^) has turned into Eq. (|l^) . Besides its nonlinearity corresponding to the 
gap the latter has a similar structure as the linear Eq. (0). One also realizes 
that for the zero density limit N/V Eq. (|lj) turns into the Schrodinger 
equation for the deuteron. This is demonstrated in Fig. ^ where we show the 
Fourier transform of tZ/p, i.e. the Cooper pair wave function for decreasing den- 
sities (Fig. ^ is taken from jl^). One clearly sees that the np Cooper pair wave 
function approaches the free deuteron wave function in the zero density limit. 
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Figure 7. Triton binding energy per nu- 
cleon as a function of nuclear density at 
T = 10 MeV, and total momentum rel- 
ative to the medium. From left to right: 
P = 0, 1,2, 3,4,5 fm-i (l|. 



Figure 8. Mott momentum per nu- 
cleon PiAott/A for triton - solid (T = 10 
MeV) and dashed line (T = 20 MeV) and 
deuteron - dashed dotted (T = 10 MeV) 
and short dashed line (T = 20 MeV) as 
a function of nuclear density ^] . 



4 Three-body reactions and bound states 

A first three-body equation for quasiparticles embedded in a medium at zero 
temperature has been given in Here we present the equation for the wave 
function 

{E- 61-62- es)ME) = [{fih - /i/2)1^2(12) + perm.] ^^iE) (16) 
together with the normalization 

E ^IWlhh + flf2f3)^l23 = (17) 
123 

where the dual wave functions are defined by by = (/1/2/3 + /i/2/3)^i23- 
An extension to finite temperatures uses the Alt-Grassberger-Sandhas 
(AGS) formahsm to treat scattering phenomena and has been given in p^ . 
For a numerical solution of the scattering problem we use a separable ansatz for 
the strong nucleon-nucleon potential 1^] . In |^ we have calculated the break up 
cross section Nd —> NNN and found considerable dependence of the deuteron 
fluctuation time on the proper treatment of the medium dependence 0, see 
also Fig. ^. The influence of this refined treatment on the relaxation times is 
treated in [|oj. The solution for the medium modified three-body bound state 
using the Yamaguchi and the Paris potential for comparison is provided in [ p^ . 
Fig. shows the triton binding energy depending on different cm. momenta as 
a function of the nuclear density and Fig. || shows the Mott momentum. Below 
the respective lines no deuterons/triton exist as bound states. 

A very interesting situation can also happen when a sharp pole at -E = 3^ 
appears in the in- medium three-body equation (p^. Then the original Fermi 
gas becomes unstable with respect to the new Fermi gas composed out of three- 
body clusters. For a very low density nuclear matter this could be a Fermi gas 
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Figure 9. Critical temperatures for the onset of quantum condensation in symmet- 
ric nuclear matter, model calculation. The critical temperature of the onset of two- 
particle pairing TJr is compared with r| for the onset of a four-particle condensate, 
as a function of the chemical potential (left) and as a function of the uncorrelated 
density (right) ^. 



of ■^He or tritons. Though this case may be somewhat academic because of the 
presence of other type of clusters, e.g. a-particles, one may think of another 
situation where this consideration may become very relevant; suppose one cools 
down a quark gluon plasma to the point that nucleons (and other hadrons) 
appear. Then the quark Fermi gas transforms into a nucleon Fermi gas. A 
relativistic equation describing this transition should have a similar structure 
as Eq. (|l^). From such an equation the critical temperature for the transition 
from quark gas to nucleon gas could be studied. 



5 Four-nucleon correlations and quartetting 

The four-nuclcon correlation is believed to play a significant role for lower 
densities and temperatures because of the strong binding energy of the a- 
particle. Exploratory calculations using a simple variational ansatz for the "^He 
wave function predict an a condensate/quartetting on top of the deuteron 
condensate/triplet pairing that leads to superfluidity 

Bose systems behave quite differently with respect to the occurrence of 
bound states. The Bose functions enhance the effective residual interaction 
that might lead to an "opposite Mott effect" , i.e. existence of bound states and 
also pairing (condensate) even if no bound state exist for the isolated case. An 
example is provided by a pion gas, were a pionic condensate may occur ||l^ 
that are discussed, e.g. in the context of neutron stars [§. 

In the framework utilized here the four body Schrodinger equation for the 
bound state problem reads 



- £1 - £2 - £3 - SiMiE) - [(/1/2 - fif2)V2{12) + perm.] ^:({E) (18) 
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with the normalization [ [TgI 

E '^r234(/l/2/3/4 - /l/2/3/4)'^r234 = , (19) 
1234 

and the dual wave functions 3'f234 = (/1/2/3/4 - /i/2/3/4)^!'234- 

The solution for if4(4/i) provides the critical temperature T4 for a conden- 
sates as a function of the chemical potential shown in Fig. ^ (left), or of the 
density as shown in Fig. ^ (right). Using a variational ansatz to model the 
four-body wave function the transition to quartetting outreaches the transition 
to isospin singlet pairing if the density is smaller than 0.03 fm~^, see Fig. ^. 
This may again be relevant in expanding nuclear matter or even in the far tail 
of the density of a nucleus. 

One may suspect that the asymmetry of nuclear matter hinders the forma- 
tion of np pairing or of quartetting. This is because extra neutrons Pauli-block 
the formation of Cooper pairs. However, once the density is so low that the 
corresponding chemical potentials turn negative, then the blocking effect, i.e. 
the Pauli principle is less effective and formation of deuterons and a parti- 
cles should not be much unfavored in asymmetric matter with respect to the 
symmetric case. 

6 Conclusion 

We have studied the formation of clusters and their condensation in Fermi 
systems using the Dyson equation approach to many-body Green's functions. 
They only depend on the cm. energy of the cluster and therefore the resulting 
equations are of the Schrodinger type. To lowest order the in-medium n-body 
equations are - compared to the isolated case - only modified by RPA phase 
space factors on the two-body matrix elements (leading to self-energy correc- 
tions and modified interactions). Such equations have early been written down 
by the Japanese school . We here investigate them with respect to the Mott 
effect and Bose-Einstein condensation in homogeneous matter. It is found that 
deuteron Cooper pairs survive longer than a-particle quartets at high densities 
but that, not unexpectedly, a-particlcs win in the low density limit. We give the 
Mott density for deuterons and tritons and discuss the Pauli principle aspect 
of the phase transition of a quark Fermi gas to a nucleon Fermi gas. Life time 
effects of the clusters are also considered. To second order Born approximation 
they behave like Fji = a {E — nep)'^ + bT^, where n is the nucleon number in 
the cluster. Therefore n-body clusters at — nep and T = have infinite life 
time. 

Although we have taken examples from nuclear physics, the exact treat- 
ment of few-body systems embedded in the medium is not restricted to nuclear 
physics. Further potential applications are possible, e.g. in the field of semi- 
conductors to treat the formation of excitons or the trion bound state in the 
dense low dimensional electron plasmas. For recent experiments see e.g. Q. 
For highly ionized dense plasmas the impulse approximation may fail for the 
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three-particle break-up cross section needed to calculate the ionization and 
recombination rates. 
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